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ABSTRACT 

The  possibility  Is  examined  of  studying  the  behavior  of  elastic 
solids  on  the  basis  of  a model  which  considers  them  to  be  made 
up  by  a very  large  number  of  elementary  particles  arranged  In  a 
regular  lattice  array.  Some  of  the  simplifications  which  are 
needed  (e.g.,  the  neglect  of  surface  tension)  to  make  the  anal- 
ysis practical  are  examined,  first  In  general  terms  and  then 
for  the  simple  case  of  a bar  In  tension. 


*The  work  described  here  was  supported  by  the  Office  of  Naval  Research 
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1,  Introduction 


The  analysis  of  the  stresses  and  deformations  of  solids  is  normally  carried 
out  on  the  basis  of  a continuum  model,  whether  classical  or  exhibiting  some  sort 
of  microstructure,  which  represents  macros cop ically  the  actual  discrete  struc- 
ture of  the  material.  Alternatively,  a discrete  model  can  be  used,  and  there  are 
then  two  basic  ways  in  which  one  might  proceed.  In  the  first,  we  start  from  the 
continuum  model,  and  replace  it  by  a discrete  one  by  means  of  suitable  partition- 
ing into  appropriate  units.  This  is  the  basis  of  very  powerful  methods  of  anal- 
ysis, such  as  those  employing  finite  elements  or  finite  differences.  The  great 
value  of  these  approaches  is  well  known  and  need  not  be  discussed  here. 

A second  approach  consists  in  starting  with  the  discontinuous  structure  and 
attempting  to  perform  calculations  immediately  on  that  basis.  This  has  been 
tried  (e.g. , [1,2]) ; in  these  works  the  body  is  considered  to  be  composed  of  a 
relatively  small  number  of  material  points,  interactions  between  them  being  con- 
sidered numerically.  To  represent  more  realistically  the  actual  material,  how- 
ever, one  would  think  that  a very  large  number  of  material  points  should  be  used — 
so  large,  in  fact,  that  a continuum  may  in  some  sense  be  approached.  It  is  un- 
likely that  a purely  numerical  approach  is  useful  in  this  connection,  because  of 
the  enormous  number  of  points  being  considered;  it  is  however  possible  that 
analytical  techniques  may  be  introduced  to  assist  in  the  calulatlons.  An  exam- 
ination of  the  latter  possibility  forms  the  subject  of  this  paper. 

The  present  approach  starts  with  the  consideration  of  a large  number  of 
material  points,  representing  atoms  and  thus  being  assumed  to  be  separated  by 
distances  of  atomic  dimensions.  Atomic  forces  are  assumed  to  act  between  them 
in  accordance  with  the  appropriate  laws  of  attraction  and  repulsion,  necessitating 
the  establishment  of  equilibrium  equations  at  each  point.  These  equations  are 
then  simplified  by  neglecting  the  effect  of  surface  tension  [3],  and  by  the  use 
of  Saint  Venant's  principle;  this  allows  a mathematical  shift  to  a continuum 
and  thus  solution  of  the  problem. 


I 
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This  approach  therefore  Involves  postponing  as  long  as  possible  the  tran- 
sition from  the  lattice  to  the  continuum  model.  This  has  the  advantage  that  the 
transition  from  one  to  the  other  can  be  examined  from  the  standpoint  of  a purely 
mathematical  approximation,  and  thus,  at  least  In  principle,  any  errors  Incurred 
can  be  easily  ascertained.  The  problem  chosen  here  Is  the  simplest  one  possible, 
namely  that  of  a bar  in  tension,  and  one  may  have  serious  doubts  whether  Che 
present  approach  can  be  extended  to  any  but  Che  most  elementary  cases.  It  Is 
nevertheless  hoped  that  the  present  analysis  may  be  of  some  Interest  In  Illus- 
trating Che  concepts  Involved  and  the  difficulties  that  might  be  encountered. 


2.  Basic  Considerations 


Consider  two  particles,  located  respectively  at  Pj^(x^,yj^,Zj^)  and 


Pj(Xj,yj,Zj)  and  separated  by  a distance 


= Let  the  force  between 


them  be  along  r^^,  l.e.,  so  chat  its  components  are 


,yry* 


(1) 


where  |F^j|  > F^j.  During  deformation,  the  particles  undergo  displacements 
and  Uj  respectively.  The  consequent  change  In  Che  components  (1)  Is  easily 
calculated,  and  gives  the  additional  force  between  the  particles,  on  the  assump- 
tion of  Infinitesimal  deformations,  as 


ij  2 r^^  \ j ij 


(2) 


r.  (u,-uj 

r“ 

ij 

It  will  be  next  assumed,  as  Is  usual  [4,5],  Chat  the  force  between  particles 


represents  Che  difference  between  an  attractive  (positive)  and  a repulsive  term,  or 


_ _ aV  _ mA nB 

j " ar, , " m+1  n+1 

‘J  'ij  'ij 


(3) 


r 


corresponding  to  a potential 


n > m > o 


(3a) 


where  A and  B are  constants.  The  term  in  parenthesis  is  eq.(2)  is  then 


n(n  -I-  2)B  _ m(m  + 2)A 


n+2 

•ij 


ia4'2 

•ij 


(4) 


It  is  now  easy,  in  principle,  to  see  how  one  should  proceed:  the  resultant 
force  on  each  particle,  exerted  by  all  other  particles,  is  set  equal  to  zero 
if  the  particle  is  Internal,  and  equal  to  the  appropriate  applied  force  if  the 
particle  is  on  the  surface  (on  the  assumption  that  only  surface  forces  are 
applied).  Practically,  however,  this  would  give  rise  to  a prohibitive  number 
of  equations  and  unknowns,  and  it  is  therefore  necessary  to  simplify  the  problem. 
To  that  end,  two  surface  effects  will  be  neglected.  The  first  of  these  is  the 
one  conventionally  disregarded  in  continuum  solutions  on  the  basis  of  Saint 
Venant's  principle,  which  allows  the  details  of  distribution  of  the  forces 
applied  over  a small  surface  are  to  be  overlooked,  since  they  practically  do  not 
affect  the  solution  outside  a thin  surface  layer. 

The  second  neglected  effect  does  not  arise  In  classical  elasticity, 
although  It  arises  in  continuum  theories  exhibiting  microstructure  [3].  It 
refers  to  the  presence  of  surface  tension,  and  Is  also  confined  to  a thin  layer 
adjacent  to  the  surface  of  the  body.  The  thickness  of  this  layer  however,  of  the 
order  of  the  atomic  dimensions,  and  therefore  much  smaller  than  that  of  the 
Saint  Venant  layer,  the  latter  being  of  the  order  of  the  linear  dimensions  of 
the  loaded  area.  Hence  surface  tension  affects  an  extremely  small  portion  of  the 
body  and  can  safely  be  neglected  if  (as  is  presently  desired)  a solution  anal- 
ogous to  the  classical  one  Is  sought. 

To  clarify  the  above  general  remarks,  the  special  example  of  a uniform  bar 


In  tension  will  be  considered. 


3,  A Bar  in  Tension 

Consider  s uniform  bar,  extending  in  the  x-dlrectlon,  and  pulled  in  that 
direction  by  an  end-force  P,  all  other  surfaces  being  free  of  applied  loads. 
Consider  first  the  forces  In  the  y-dlrectlon  on  a generic  1-Ch  particle; 
equilibrium  requires  that 


^ 6F, - o 


Ijy 


for  all  1 


(5) 


all  J 


where  5P^jy  is  the  y- component  of  Assume  for  slidpllclty  of  discussion 


that  the  particles  occupy  the  vertices  of  a simple  cubic  lattice,  and  that 
y^constant  represents  one  of  the  free  faces  of  the  bar.  Consider  first  the 
portion  of  the  summation  which  refers  to  the  y-direction:  as  one  proceeds 
further  and  further  from  a face  into  the  Interior  of  the  body,  more  and  more 
planes  of  particles  appear  In  the  summation.  This  gives  rise  to  the  surface 
effect  previously  referred  to,  and  which  is  limited  to  very  few  such  planes;  to 
disregard  It  Is  equivalent  to  neglecting  the  presence  of  a physical  surface, 
except  as  a place  where  loads  may  be  applied.  In  other  words  the  body  is  math- 
ematically considered  to  extend  beyond  Its  actual  surface,  sufficiently  far 
that  the  same  type  of  equilibrium  equations  hold  at  all  planes  regardless  of 
their  distance  from  the  surface.  The  same  argument  holds  for  all  surface  of 
the  body,  and  eq.(5)  for  our  problem  can  then  be  written  as: 

R R R 

I I I 


ijy 


for  all  1 


(6) 


yj- 


Zj-  -R 


where  R Is  a distance  chosen  to  be  sufficiently  large  to  Include  the  entire  bar 
and  what  additional  portion  Is  needed  to  eliminate  the  surface  effects;  it  will 
be  seen  that,  for  ease  of  calculations.  It  may  at  times  be  taken  as  infinite. 


J 
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For  any  value  of  1,  the  y-portlon  of  the  summation  can  now  be  separated  into  the 
two  regions  y^  S y^  and  y^  ^ ^1*  since  no  forces  In  this  direction  are  any- 
where applied,  each  of  these  regions  must  separately  yield  a zero  resultant,  or 
(with  symmetry  In  the  other  directions): 


R ^1  R 

I I "lly-” 

Xj-o  Yj-  -R  Zj-o 


‘I 


I 


I 

Z ,*o 


o 


(7) 


The  treatment  of  equilibrium  In  the  z-dlrectlon  Is  entirely  analogous;  so  Is 
that  for  the  x-dlrection  except  for  the  Introduction  of  the  applied  force  P,  or 


R 


‘I  I 


R 
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ijx 


xj-  -R  y-o  Zj^ 


(8) 


R 


‘Z  I 


R 

r* 


L fi^ijx 


-p 


In  this  manner  the  local  surface  effects  have  been  eliminated,  but  the  presence 
of  applied  surface  loads  has  been  accounted  for. 

Assume  now, either  by  similarity  to  the  classical  solution  or  as  an  inde- 
pendent semi-inverse  assximptlon,  and  on  the  basis  of  Saint  Venant's  principle, 
that  the  displacement  components  of  a generic  point  P^  can  be  taken  to  be 


"i  “ “i  • 


^4 


Wi  - bz^ 


(9) 


where  a and  b are  constants.  The  components  of  then 


6F 


ijx 


b^^- 

drij 


(a  - b)(- 


dr 


ij 


F ^ 
'ij''  r^ 


(10) 


where 


5ij  - - "i  • ^ij  “ • qj 


(11) 


The  second  ones  of  eqs.  (7)  and  (8)  then  become,  respectively, 
bSj^  + (a  - b)S2  » o 


(12) 


as 


^ + 2(b  - a)S2  - P 


where 


'.-I  I Ift-I  I 

§“o  r^h  ^=o  5=h  TT=o  ^=o  ^ 


^ ^ ^ SF  P 2 

1-0  irh  C=o  ^ ^ ij 


(13a) 


(13b) 


and  where  h - ~ Is  the  side  of  the  basic  cube  in  the  lattice.  The 

other  forms  of  S2»  obtained  by  permutation  of  5,ti  and  have  not  been  shown  in 
eq.  (13b)  for  the  sake  of  brevity,  and  the  subscripts  ij  have  been  dropped  from 
these  quantities  for  simplicity  of  writing.  It  then  appears  that  the  first  of 
eqs.  (7)  and  (8)  are  also  satisfied,  because  use  of  (12)  automatically  insures 
that  the  equilibrium  equations  are  satisfied  for  all  values  of  i. 

Solution  of  (12)  gives 


P(S2  - S^) 


PS, 


Si(3S2  - S^) 


• S^(3S2  - Sp 

In  particular,  we  may  define  Poisson's  ratio  as 
. S„ 


(14) 


a S,  -S2 


(15) 


There  now  remains  the  evaluation  of  and  S^:  this  may  be  performed  either 
in  a direct  numerical  fashion,  or  in  an  approximate  analytical  manner  by 
replacing  the  sum  over  discrete  points  by  integrals  over  a continuum.  The  latter 
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approach  will  now  be  followed;  it  exemplifies  the  previously  discussed  trans- 
ition to  a continuum  model  at  a late  stage  of  the  analysis. 


4.  Evaluation  of  Sums 

The  trapeziodal  rule,  usually  employed  to  approximate  integrals  by  means  of 
sums,  may  be  conversely  used  to  approximate  and  present  sums  by  means  of  integrals 
Then,  provided  R is  large  enough. 


'i  ~ f f 5 TT  r d r d 5 

1 J J dr  ^ o o ^ 

u u 


^ ^ ^ r^F  F 2 

j j 

h h h ^ ^ ij 

2 2 2 

where  r = ri  + C • Then,  with  (3)  and  (4),  we  have 
o 


g-i  R)“‘^  Cv^h)' 


>(Hi) 


(JZR)^'^  Cj2h)' 


f,  (°iR^+2hS 

(rW)  ^ (m-2)  2(m-2)(v^h) 


(nR^+2h^')  _ (n-h2) 

2(n-2)(v^  h)""^ 


Evaluation  of  these  quantities  requires  a knowledge  of  m,n,A,B,  and  h,  as 
well  as  a choice  of  R.  We  can  however  derive  some  general  information,  partic- 
ularly about  Poisson's  ratio,  without  such  detailed  knowledge,  but  by  inserting 
eqs.(17)  into  (15),  and  then  simplifying  the  resulting  expression  on  the  basis  of 
the  inequalities 

n-m 


R ‘ ^ ’ B 


1 Ah“  , 
<*  1 ; — — < 1 

BR 


(18) 


and  recalling  that  m < n.  Three  cases  must  now  be  distinguished,  namely: 

(m  > 2,  n > 2),  (m  < 2,  n > 2),  and  (m  < n < 2);  then  (18)  permits  the  expres- 
sions obtained  for  Poisson's  ratio  from  eq.(15)  to  be  reduced  to: 


\) 


n -H  2 


3(n  -f-  2) 


for  n > 2,  any  m < n 


(19a) 


m 


m + 2 


for  m < n < 2 


(19b) 


It  is  easy  to  see  that,  from  eq.  (19a),  (1/3)  < v < (1/2),  while  from  eq.  (19b) 
we  have  o < v < (1/2),  Thus,  in  agreement  with  experimental  observation,  v never 


exceeds  1/2,  but  any  positive  value  less  than  1/2  is  possible.  Most  metals,  for 


example,  fall  in  the  former  category  [5],  Analyses  such  as  the  present  one  can 

4 

be  used  to  determine  m or  n from  a knowledge  of  Poisson's  ratio,  or,  conversely, 
to  predict  v on  the  basis  of  a postulated  atomic  behavior.  The  adaptation  of  the 
concepts  discussed  here  to  more  complicated  problems  will  of  course  have  to  await 
further  work  along  these  lines. 
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